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PREFACE 

!nals  work  is  primarily  concerned  with  the  interaction  between 
the  flow  of  an  ionized  gas  and  a  magnetic  field.  The  particular 
problem  treated  corresponds  very  closely  to  conditions  existing  in 
exploding  wire  experiments,  and  should  prove  useful  in  interpreting 
the  results  of  such  experiments. 
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SUMM/iJ?Y 


A  similarity  solution  is  obtained  for  the  flow  behind  a  very 
strong  (in  the  hydrodynamic  sense)  cylindrical  magnetohydrodynamlc  shock 
wave  produced  by  the  sudden  release  of  energy  along  a  line  of  infinite 
extent  in  a  plasnsa^  -  The  plasma  is  assvuned  to  be  an  ideal  gas  with  in¬ 
finite  electrical  conductivity,  and  to  be  permeated  by  the  azimuthal 
magnetic  field  of  a  line  current.  It  is  shown  that  it  is  of  critical 
importance  to  take  into  account  the  ambient  magnetic  pressure,  no  matter 
how  small.  It  is  found  that,  to  preserve  similarity,  the  external  circuit 
is  required  to  maintain  a  constant  axial  current;  this  result  also  appears 
in  the  related  problem,  treated  by  Greenspan,  where  the  ambient  plasma 
la  nonconducting.  It  is  shown  that  this  boundary  condition  on  the  circuit 
Implies  an  exchange  of  energy  at  t  =  0  between  the  external  circuit  and 
the  plasma.  When  this  energy  is  taken  into  account,  the  dependence  of  the 
shock  speed  on  the  explosive  energy  can  be  obtained.  This  dependence  is 
determined  as  a  function  of  the  ambient  magnetic  field  both  for  the  present 
case  and  for  Greenspan's  case,  and  Inteii'stlng  differences  are  noted.  The 
fraction  of  explosive  energy  which  appears  as  mechanical  energy  is  also 
calculated  for  the  two  cases,  and  again  significant  differences  are  found 
to  exist.  Other  differences,  with  possible  experimental  consequences,  are 
also  discussed.  The  distributions  of  flow  quantities  in  the  neighborhood 
of  the  axis  In  the  present  case  la  found  to  be  different  from  both  the 
ordinary  blast  wave  and  Greenspan's  case;  the  most  pronounced  differences 
arc  In  the  density  distribution  for  large  magnetic  field  and  in  the  pressure 
distribution  for  ani*  magnetic  field  whatever. 
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I .  IWTRODUCTICW 


It  l8  veil  knovn  that  a  siinllarlty  solution  exists  for  the  flow 
behind  a  very  strong  cylindrical  shock  wave  produced  by  the  sudden  re¬ 
lease  of  a  finite  amount  of  eneigy  per  unit  length,  along  a  straight 
line  of  Infinite  extent  In  a  uniform  gas.^^^  A  corresponding  magneto- 
hydrodynamic  problem  also  has  the  same  similarity.  '  Here  the 

shocked  gas  Is  an  electrically  conducting  plasma  which  Interacts  with  an 
applied  magnetic  field.  The  gas  Is  assumed  to  be  Initially  In  a  purely 
azimuthal  magnetic  field  produced  by  a  constant  line  current,  the  return 
for  which  Is  a  concentric  cylindrical  conductor  of  very  large  radius.  The 
flow  Is  again  produced  by  the  sudden  release  of  energy  along  the  axis  (by 
exploding  the  wire  carrying  the  current,  for  example). 

The  flows  considered  here  are  those  for  which  the  conductivity  of 

(U) 

the  shocked  plasma  Is  infinite.  For  this  class  of  flows,  as  Greenspan^  ' 
has  pointed  out,  two  limiting  cases  arise.  In  the  first  case,  the  quiescent 
gas  has  zero  conductivity  but  the  gas  behind  the  shock  has  Infinite  con¬ 
ductivity.  The  magnetic  field  across  the  shock  Is  continuous  and  the 
boundary  conditions  at  the  shock  are  the  oi^dlnary  hydrodynamic  Jump  condi¬ 
tions  for  a  strong  shock.  This  situation  arises  If  the  magnetic  shock 
transition  zone  Is  much  wider  than  the  viscous  transition  zone,  which  Is 
replaced  by  a  shock  Jump.  The  entire  effect  of  the  magnetic  field,  In 
this  case,  la  contained  In  Its  Interaction  with  the  flow  behind  the  shock. 
This  limiting  case,  hereinafter  referred  to  as  "modified  hydrodynamic," 
has  been  discussed  by  Greenspan.  In  the  second  limiting  case,  the  magnetic 
field  Is  discontinuous  across  the  shock,  and  the  boundary  conditions  at 
the  shock  are  those  appropriate  to  a  magnetohydrodynamlc  shock.  This 
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sltuation  results  if  both  the  magnetic  ana  viscous  transition  zones  are  , 
very  thin  and  if  the  quiescent  gas  is  also  a  perfect  conductor  (as  a  result 
of  pre-ionization,  for  example).  In  this  case,  hereinafter  referred  to  as 
"pure  MHD, "  the  magnetic  field  affects  conditions  at  the  shock  itself  as 
veil  as  those  behind  it. 

In  this  paper,  the  pure  MHD  flow  is  considered  in  detail.  This  case 
was  first  discussed  by  Pal,^^^  but  the  results  obtained  below  are  found  to 
disagree  significantly  with  his.  The  discrepancy  is  shown  to  arise  from 
an  approximation  made  by  Pai  to  the  pressure  Jump  condition  at  the  shock. 
The  partic\ilar  approximation  was  to  neglect  the  magnetic  pressure  in  front 
of  the  shock.  While  it  is  true  that,  fur  small  values  of  this  magnetic 
pressure,  the  errors  introduced  at  the  shock  are  small,  nevertheless  the 
inclusion  of  such  n  pressure,  no  matter  how  smill,  produces  a  profound 
effect  at  the  axis  and,  concomitantly,  in  the  energy  content  of  the  gas. 

The  question  of  the  energy  content  of  the  gas  is  also  considered  below, 
for  the  modified  hydrodynamic  as  well  as  for  the  pure  MHD  case.  It  is 
Shown  that  in  both  cases,  in  contrast  to  the  ordinary  blast  wave  solution, 
the  Increment  in  the  total  (gas  uyna.-^ic  plus  raignetic)  energy  of  the  gas 
is  not  equal  to  the  energy  suppliea  by  tne  explosion.  The  two  cases  differ 
markedly,  however,  in  that  in  the  pure  MHD  case  the  increment  is  infinite, 
whereas  in  tne  modified  hydrodynamic  case  it  is  finite  but  lens  than  the 
energy  supplied  by  the  explosion.  The  difference  between  the  Increment 
ana  the  explosion  energy  represents  n  concent mt ion  of  magnetic  energy, 

In  one  case  supplied  by,  ana  in  the  other  cane  removed  by  the  external 
circuit  In  aalntaining  the  constant  current  required  by  the  similarity 
solutlonr- — TE"il  "possible  in  both  eases  to  calculate  the  contribution  from 
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the  external  circuit,  and  thereby  to  relate  the  shock  speed  and  the 
mechanical  energy  in  the  gas  to  the  energy  supplied  by  the  explosion. 

The  analysis  reveals  an  Interesting  difference  between  the  two  cases. 

In  both  cases  only  part  of  the  energy  supplied  by  the  explosion  appears 
as  mechanical  energy,  the  remainder  being  converted  to  magnetic  energy. 

In  the  pure  MHD  case,  this  converted  part  of  the  explosion  energy  appears 
in  the  compressed  magnetic  field  behind  the  shock,  while  in  the  modified 
hydrodynamic  case,  it  is  delivered  Insteaa  to  the  external  circuit. 

It  seems  reasonable  to  believe  that,  despite  the  divergence  mentioned 
above,  the  results  obtained  have  physical  significance.  The  divergence 
appears  because  of  the  unphysical  assumption  of  a  dimensionless  axial  con¬ 
ductor.  One  would  expect  that  in  any  actual  experiment,  where  the  axial 
conductor  has  a  finite  radius,  the  true  flow  wouIq  approach  that  given 
by  the  similarity  solution  at  alstances  from  the  axis  much  greater  thaji 
the  radius  of  the  conductor.  One  wouIq  also  expect  that  any  relations  ob¬ 
tained  between  finite  quantities,  such  as  the  dependence  of  the  shock  speed 
on  the  explosion  energy  or  the  ratio  of  mechanical  energy  to  explosion  energy, 
WouIq  also  be  approximately  valid.  The’''*  are  questions  which,  of  course, 
must  ultimately  be  decided  experimentally. 


II .  FORMULATION  OF  THE  PROBLEM 


A  solution  of  the  magnetohydrodynaalc  equations  is  desired  for  the 
flow  behind  a  cylindrical  shock  vave  produced  by  the  sudden  release  of  a 
finite  amount  of  energy  per  unit  length,  along  an  axis  of  infinite 
extent  in  a  plasma.  The  plasma  is  assumed  to  be  an  ideal  gas  with  infinite 
electrical  conductivity;  viscosity  and  heat  conduction  are  neglected.  An 
azimuthal  magnetic  field  is  assvused  to  exist  initially  in  the  plasma.  As 
a  consequence  of  the  assumption  of  infinite  conductivity,  the  magnetic 
field  in  the  plasma  remains  azimuthal.  Under  the  assumptions  made,  all 
quantities  are  functions  only  of  the  radius  r  and  time  t,  and  the  motion 
of  the  gas  is  governed  by  the  following  equations: 


r|f.^(«u) 


St  dr  p  dr  p  dr 


fo 


S. 

Pur 


sr 


0 


0 


(continuity) 

(momentum) 

( induction) 
(entropy) . 


(la) 

(lb) 


(lc) 

(ld) 


In  these  equations,  p  la  the  gas  pressure,  p  is  tue  density,  u  the  radial 
velocity,  and  Bg  the  azimuthal  magnetic  field.  The  specific  heat  ratio  y 
is  assumed  constant  and  the  magnetic  permeability  u  is  taken  to  be  that 
of  free  space  (MKS  units  are  used  throvighout ) .  The  gas  is  assumed  to  be 
initially  in  a  uniform  state  p^,  p^.  The  initial  azimuthal  magnetic  field 
is  assumed  to  be  that  produced  by  a  constant  current  I^,  flowing  along  the 
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Entering  into  the  initial  and  boundary  conditions  of  the  problem  are 

certain  dimensional  constants;  these  are  the  initial  pressure  p^,  the 

Initial  density  the  initial  axial  cxxrrent  and  the  energy  11b- 

2 

erated  per  unit  length.  Since  the  dimensions  of  the  quantity  are  the 
same  as  those  of  only  three  of  the  four  constants  are  dimensionally 
independent.  If,  in  addition,  the  liberated  energy  Is  large  enough  to 
produce  a  veiy  strong  shock,  the  initial  pressure  p^  may  be  neglected,  and 
only  two  constants  with  Independent  dimensions  enter  into  the  problem— p^ 
and  E^.  As  in  the  ordinary  blast  wave  case,  the  only  combination  of  length 
and  time  that  can  be  formed  from  these  constants  is  of  the  form 


time 

It  now  follows  from  general  considerations  that  the  system  of  equations, 
(1),  possesses^ Similarity  solution;  this  solution  is  of  the  form 


P  ■  0^  (q) 

(2a) 

U  -  1  “2^  (q) 

(2b) 

p  ■ 

(2c) 

1 

'O  ■  t 

(Pd) 

where  ^ ^ and  S  are  non-dimensional  functions  of  a  non-dimensional 
variable.  The  only  non-dimensional  variable  which  can  be  formed  from  the 
variables  r  and  t  and  the  dimensional  constants  of  the  problem  is  the 


variable 
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T  -  ~ 
c  t 


vhez«  c  is  some  characteMstlc  parameter  with  the  dimensions  of  (length)  / 


Zt  Is  convenient  to  Introduce  as  independent  variables,  Instead  of  r 


and  t,  the  stream  function  ^  defined  by 


dr  p  dt  0 

O  0 


and  the  variable  defined  as 


r 


Physically,  the  value  of  corresponding  to  any  (r,t)  is  the  ratio  of  the 
mass  between  the  axis  end  the  radius  r  at  tlmettothe  mass  Initially  in 
the  same  volume.  The  value  of  0  at  the  shock,  therefore,  is  0^  ■  1.  A 
line  of  constant  0  is  also  a  line  of  constant  n  in  the  (r,t)  plane. 

Zn  terms  of  the  variables  (t»0),  'vds.  (?)  take  the  form 


P  ■  Oq  c(0) 


u  ■  c 


P  . 


The  quantity  q  defined  by  Eq.  (5)  Is  proportional  to  the  square  of 
that  defined  by  Greenspan. 
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vhere  the  quantities  a,  U,  P,  and  3  are  non-diaensional  functions  of  the 
single  (non-dimensional)  variable  f).  The  parameter  c  in  Eq.  (6)  is 
related  to  the  speed  of  the  shock;  in  general, 

c  -  A  ^  (7) 


where  c  is  the  shock  speed  and  A  is  an  arbitrary  constant*  In  the  present 

case,  it  is  convenient  to  choose  A  »  1. 

It  is  necessary  for  the  initial  conditions  to  be  compatible  with  the 

assumed  form,  f6),  of  the  solution.  The  plasma  is  assumed  to  be  initially 

at  rest,  to  have  a  uniform  density  p  ,  a  uniform  pressure  p  ,  and  to  be 

o  o 

permeated  by  the  azimuthal  magnetic  field  of  a  constant  line  current  I  . 

0 

The  first  two  conditions  are  obviously  expressible  in  the  desired  form; 
viz., 

-  0  (8a) 

o 

and 


o  ■  1 
o 


(8b) 


respectively,  while  the  last  condition  ca;'*  be  written  in  the  form 

1/? 

I  uC  1 


p  c 


(8c) 


For  the  condition  of  uniform  pressure  to  be  compatible  with  (6),  however, 
It  Is  necessary  to  assume  that  the  initial  pressure  Is  negligible;  i.e. 


P 


0 


0  . 


(8d) 


This  Is  the  usual  blast  wave  assumption  of  a  very  strong  shock. 
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0  (mom^ntuffl)  (l*?b) 
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physical  variables,  these  conditions  may  be  written: 


P  u 

_o  _  _ s 

P-  “  ’  c 


p.  *  h 


B  u 

-£  .  1  - 

B.  C 

8 


0  u  c 
0  s 


(c^  -  <2|i  cUg)(c-Ug)^ 


B" 

0 

PP. 


(=^ 


^  cu  +  ■?  u  ) 
2  s  2  s'^ 


( induction) 

(12c) 

(entropy)  . 

(12d) 

s  the  boundary  condl- 

shock.  In  terms  of 

( continuity) 

(13a) 

(momentum) 

(l^b) 

( induction) 

(IJc) 

(wave  speed) 

(13i) 

where  the  subscript  "s"  denotes  quantities  immediately  behind  the  shock 

and  B  is  the  value  of  B-  immediately  ahead  of  the  shock.  The  quantity 
c  ^ 

c  is  again  the  shock  speed,  given  by  (7'  with  A  «  1.  The  wave  speed  equa¬ 
tion,  (13d),  is  derived  by  combining  the  usual  energy  equation  with  the 
continuity,  momentum  and  induction  equations. 

It  should  be  noted  that  the  magnetic  pressure  ahead  of  the  shock  has 
teen  retained  in  the  momentum  equation,  (I3b).  AJthough  it  is  necessary 
to  neglect  the  ambient  gas  pressure  in  this  equation  to  preserve  similarity, 
the  same  is  not  true  of  the  ambient  magnetic  pressure  which  does  have  the 
proper  form.  This  term,  which  was  neglected  by  I'al,  can  therefore  be 
retained  ana,  as  will  be  shown  below,  should  be  retalnea. 
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In  teras  of  the  non-dimensional  variables  defined  by  Eq.  (6),  the 
shock  relations  become: 


^  -  1  -  U 

C  8 

(14a) 

s 

p.  *  1  - 1)  -  U3 

(I4b) 

B. 

P.  ’  ^ 

(I'^c) 

(1  -  U,)  -  (1  -|u,)  • 

(l^d) 

The  quantity  is  given  by  Eq.  (8c),  and  is,  physically,  the  ratio 
of  the  Alfye^'n  speed  in  the  gas  immediately  ahead  of  the  shock  to  the 
shock  speed  (Alfven  number). 

If  is  considered  a  parameter,  Eqs.  (14)  define  a  one-parameter 
family  of  shock  curves.  For  a  given  value  of  p^,  is  first  determined 
as  a  root  of  the  quadratic  equation,  (l4d).  One  then  obtains  o  and  p 

directly  from  Eqs.  (l4a)  and  (l4c),  respectively,  and  then  P  from  Eq. 

8 

(I4b). 

Three  algebraic  integrals  of  the  four  equations  of  motion  may  be 
obtained.  The  entropy  equation,  (12d),  hss  the  immediate  integral 

p 

•  a  •  (constant)  (15) 

from  which  it  follows  that  the  entropy  is  constant  along  a  streamline.  A 
second  Integral  owy  be  obtained  by  combining  Eqs.  (IPa)  and  (12c):  namely 

•  2b  ■  (constant)  (l6) 

o 

which  expresses  the  fact  that  the  field  is  "froien"  into  the  fluld--a 
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well-known  consequence  of  the  assumption  of  infinite  conductivity.  A 
final  algebraic  integral,  obtained  by  combining  Eqs.  (I2a)  and  (12b),  is 

a  ^  i  ^  It  ^  (^7) 

This  so-called  "energy  Integral"  results  from  the  fact  that,  with  the 
assumed  similarity,  the  total  energy  is  constant  not  only  for  the  entire 
flow  field,  but  also  between  any  two  similarity  lines. 

To  establish  the  connection  between  the  variables  used  here  and 
those  of  Oreenspan,  one  additional  relation  may  be  derived;  from  Eqs.  (ll) 

5/2 

and  (lUa),  the  value  of  q  at  the  shock  is  t)  >  2'^'  .  It  then  follows 

8 

from  Eq.  (3)  that  the  location  of  the  shock  is  at  r  (t)  ■  2^^\c*t)^^^. 

The  three  Integrals  (15)>  (16)  and  (17)  may  now  be  used  to  eliminate 
the  functions  P,  fi  and  U  from  Eqs.  (12).  The  system,  (12),  is  thereby 
reduced  to  the  single  first  order  equation 

2S  .  ^  -  (in)  [k  .  b;(S.g)] 

to  0(1-5)  .  b!(55  -  !•)  ■►  k 

where 

5 • I  •  (W) 

If  neither  a  nor  b  is  zero,  no  further  analytic  integration  appears 
possible.  In  the  general  cnee,  then,  the  problem  is  reduced  to  the 
numerical  integration  of  a  single  ordinary  equation. 
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III.  GEflERAL  REMARKS 

A  few  general  remarks  can  be  made  in  connection  with  Eq.  (l8). 

Both  cases  mentioned  In  the  Introduction,  l.e.,  pure  MHD,  and  modified 
hydrodynamic,  as  well  as  the  ordinary  "blast  wave,  appear  as  special  cases 
of  Eq.  (l8).  The  solution,  in  each  case,  is  that  integral  curve  of  Eq. 

(l8)  which  starts  at  the  point  appropriate  to  the  particular  shock  and 
terminates  at  the  point  which  corresponds  to  the  axis.  (Each  point  on 
the  integral  cxirve  maps  into  a  line  of  constant  q  in  the  (r,t)-plane.) 

In  every  case,  the  point  corresponding  to  the  shock  is  a  regular  point  of 
the  differential  equation.  Since  only  a  single  integral  curve  may  pass 
through  a  regular  point,  the  solutions  obtained  in  this  manner  are  unique. 

A  point  of  central  importance  in  the  analysis  is  the  dependence  of 
the  solution  on  conditions  at  the  shock.  The  boundary  conditions  at  the 
shock  enter  into  Eq.  (l8)  through  the  constants  a,  b,  and  k,  which  are 
determined  by  evaluating  the  integrals  (15)»  (l6),  and  (17)  at  the  shock. 
Because  of  the  non-linearity  of  Eq.  (l8),  these  constants  affect  the  entire 
flow  behind  the  shock  in  a  very  complicated  way.  The  most  striking  effect, 
however,  is  on  conditions  at  the  axis.  ,^n  axis  of  symmetry  always  corre¬ 
sponds  to  a  singular  point  of  such  a  differential  equation,  l.e.,  to  a 
point  at  which  both  the  numerator  and  denominator  vanish.  The  behavior  of 
physical  quantities  at  the  axis  is  very  sensitive  to  the  location  of  this 
singular  point,  which,  in  turn,  is  determined  by  the  relationships  among 
the  constants  a,  b,  and  k  prescribed  by  conditions  at  the  shock.  As  will  be 
seen  below,  this  sensitivity  to  the  location  of  the  singular  point  manifests 
itself  in  a  profound  difference  at  the  axis,  eapeclally  in  the  pressure 
alstrlbution,  between  the  pure  KHD  and  modified  hydrodynamic  canes. 
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IV.  OM)INA;<Y  BLAST  VAVE  SOL'JTION 


For  purposes  of  illustration  and  comparison,  the  ordinary  blast 
wave  Bolucion  in  the  present  variables  is  briefly  sketched  belov.  The 
blast  vave  solution  is  recovered  in  the  limit  of  zero  magnetic  field. 
The  appropriate  shock  conditions  are  Sqs.  (lU)  vith  =  0.  These, 
applied  to  Eqs.  (15)  through  (17),  yield  a  »  b  =  k  =  0;  the 

integrals  (15)  end  (17)  therefore  become 


yfl  '7*1' 


1  rr^ 

a  U 


*  -  O  p  =  0  , 


respectively,  with  5  defined  by  Eq.  (19)* 

The  singular  point  of  Eq.  (l8)  corresponulng  to  the  axis  in  this  case 
is  the  saddle  point  cr  =  0,  5  =  —2^  .  The  behavior  of  physical' quantities 
in  the  neighborhood  of  the  axis  cun  be  ueduced  from  un  analysis  of  Eq.  (l8) 
In  the  neighborhooG  of  tae  singular  point.  However,  In  this  case  E4.  (l8) 
ern  be  Integratea  nnalyticr.liy;  the  solution  satisfying  the  boununry  condi¬ 
tions  at  the  shock  is 


Equations  (20)-(rr)  give  the  complete  solution,  in  parametric  form, 


of  the  equations  of  motion.  The  parameter  5  decreases  monotonlcnl ly  from 

the  value  nt  the  snock  to  the  va.lue  —2—  nt  the  rixis. 

7*1  7-1 
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It  can  be  seen  from  these  equations  that  P  and  U  both  become  zero 

Lizzll 

y  2 

at  the  axis,  with  P  o^  and  U  ~  a  •  Since  0  ■  o  5  —  a  at  the  axis, 

it  follows  from  Eq.  (10a)  that  t  -^o^t,  and  consequently  from  Eqs.  (6) 
for  the  physlcaJ.- variables,  that  o  -•  0,  u  -•  0,  and  pt  -•  constant  at  the 
axis*  These  are  well-known  properties  of  the  blast  wave  solution. 
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V.  PURE  MHP  CASS 

In  the  pure  MHD  case,  the  boundarj'  conditions  at  the  shock  are 
given  by  Eqs.  (1^^)  vith  0  <  <  1.  (The  density  ratio  across  the  shock 

decreases  with  increasing  becoming  unity  when  «  1.)  These  shock 
conditions,  applied  to  Eqs.  (l6)  and  (17),  yield 

b-k-|p^.  (25) 

When  these  values  of  b  and  k  are  substituted  into  Eq.  (l8),  it  is  found 
that  a  higher  order  singular  point  a  *  0,  5-1  corresponds  to  the  axis. 

The  paniraeter  5  decreases  monotonically  from  its  value  at  the  shock  to 
its  value  at  the  axis. 

An  interesting  conclusion  can  already  be  drawn  from  the  location  of 
the  singular  point.  It  is  easy  to  show,  from  the  various  defining  relation- 
ships  and  the  integral  (l6)  with  t  *  -x  ,  that  the  magnetic  field  behind 
the  shock  is  given  by 

BQ(r,t)  -  ^  I'  (r,t)  .  '  (PU) 

The  factor  multiplying  the  parameter  {  in  Eq.  (2U)  is  clearly  the  Initial 
magnetic  field.  -ThvTs,  as  might  be  expected,  the  magnetic  field  is  every¬ 
where  compressed  between  the  shock  and  the  axis.  At  the  axis,  however, 
the  magnetic  field  maintains  its  initial  value.  This  implies  that,  to 
preserve  similarity,  the  external  circuit  must  maintain  n  constant  current 

(The  sane  result  was  obtained  by  Oreenspan  in  the  modified  hydrodynamic 
case.)  Cne  might  expect  n  priori  that  u  flow  with  the  same  similarity 


would  be  produced  if  the  current  in  the  wire  were  instantaneously  increased 
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to  a  new  constant  value  since  such  a  boundary  condition  introduces 
no  new  dimensional  constants  and  is  expressible  in  a  form  compatible  with 
the  similarity.  This,  however,  is  apparently  not  the  case;  only  the 
maintenance  of  a  constant  current  preserves  similarity. 

Because  of  the  high  order  singularity  of  Eq.  (l8)  at  a  »  0,  ?  «  1,  the  de 
pendence  of  the  other  variables  on  5  in  the  neighborhood  of  the  axis  is 
most  easily'determlned  from  the  continuity  equotlon,  (12n).  An  analysis 
of  this  equations  shows  that,  in  the  neighborhood  of  o  »  0,  I  •  I, 

where  U  is  a  constant.  It  then  foUows  from  equations  (17),  (10a),  and 
o 

(11)  that 


2 


0  ■  *  ^)^ 

(26a) 

ill 

t  (5  -  1)’"^  exp 

(26b) 

.«21_ 

n  ■  -  1)  e^P  [ 

fpiT 

(26c) 

where  ■ 
equations 


,  nnd  r.nd  are  constants.  The  last  of  these 
verifies  that  this  singular  point  indeed  corresponds  to  the 


axis  (t  ■  0). 

The  behavior  of  the  physical  variables  in  the  neighborhood  of  the 
axis  can  now  be  deduced  from  the  defining  equations,  (6),  and  the  inte¬ 


grals  (15)  and  (16).  The  results  are; 
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P  - 


-  In  Ti  ^ 


T7^ 


(27a) 


-  1  2  r  1 

“  •  (->  17^7  ’’  L‘  “■"  ''J 


-  1 


(27b) 


-  ^0°  8b  -1  r  -| 
P  •  —  Tt^  ’’  L' 


(27c) 


-  Pq®  -1 

^  _  2b  Ti  ^ 


(27d) 


Afl  in  the  ordinary  blast  wave  solution,  the  density  and  velocity / 
approach  zero  at  the  axis.  The  dependence  on  however,  is  quite  diff¬ 
erent  in  the  two  cases,  the  logarithmic  dependence  being  peculiar  to  the 
pure  MHD  case.  The  density,  in  particular,  approaches  zero  much  more 
slowly  in  the  pure  MHD  case,  especially  for  the  larger  values  of  (Fig. 

1).  The  pressure  is  markedly  different,  becoming  infinite  at  the  axis 
rather  than  asymptotically  constant  as  in  the  blast  wave  solution  (Fig.  3)* 
This  infinity  will  be  discussea  further  in  connection  with  the  energy 
content  of  the  gas. 

The  results  of  the  numerical  integration  of  Eq.  (l8)  for  various 
2 

values  of  p^,wlth  y  •  l.t,  are  shown  in  Figs.  1-1*.  it  is  clear  that  the 
presence  of  a  magnetic  field  ahead  of  the  shock  has  an  Important  effect  on 
conditions  ntthe'shock  (and  everywhere  tehin'i  it)  as  well  as  those  at  the 
axis.  The  compression  across  the  shock  decreases  from  its  strong-shock  value 
of  for  zero  magnetic  field  (p^  ■  0)  to  a  value  of  unity  for  •  1.  The 

velocity  and  gas  pressure  at  the  shock  vary  accordingly,  decreasing  from 
their  strf'ng-shock  values  to  zero  over  the  range  of  p^.  (The  similarity 
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solutlon  Is  valid,  of  coxirse,  only  so  long  as  the  shock  remains  a  strong 

shock  In  the  hydrodynamic  sense,  that  Is,  so  long  as  the  gas  pressure  ahead 

of  the  shock  Is  negligible  compared  to  that  behind  It.)  The  magnetic 

pressure  at  the  shock,  on  the  other  hand,  exhibits  a  maximum  at  a  value  of 
2 

between  0.3  and  O.U.  The  appearance  of  such  a  maximum  was  noted  In  the 

Inverse  pinch  effect.  '  The  existence  of  this  maximum  depends  only  on  the 

Jump  conditions  at  the  shock,  which  are  the  same  as  in  the  inVerse  pinch. 

The  Importance  of  Including  the  ambient  magnetic  pressure  In  the 

shock  relations  Is  illustrated  by  comparison  of  the  above  results  with 

those  of  Paij-where" this  term  was  neglected.  In  the  first  place.  If  the 

ambient  magnetic  pressure  Is  neglected  In  Eq.  (l^b),  the  Alfv^n  number, 

0^,  Is  allowed  to  become  arbitrarily  large.  This  Is  In  contradiction  to 

the  limits  0  <  B  <1  obtained  above  for  a  pure  MHD  shock.  In  the  second 
*  o 

place,  all  hydrodynamic  quantities  at  the  shock  become  Independent  of  0  , 

o 

always  attaining  their  strong-shock  values.  It  can  be  seen  from  Figs.  1-3 

that  this  is  a  reasonable  approximation  only  for  very  small  values  of  B  . 

o 

The  same  error  is  Introduced  Into  the  magnetic  field  ratio  across  the 
shock  whichr^tlT Pal's  conditions,  is  always  -2^  but  which,  with  Eqs.  (lU), 
decreases  from  this  value  to  unity  at  0^  «  1.  Finally,  and  what  Is 
responsible  for  an  Important  discrepancy  even  for  the  smallest  values  of 
0^1  Pal's  conditions  lead  to  different  values  for  the  constants  a,  b, 
and  k  of  Eqs.  (13)-(17)i  and  thereby  alter  the  location  of  the  singular 
point  of  Eq.  (18)  which  corresponds  to  the  axis.  With  Pal's  conditions, 
the  singular  point  is  located  at  a  «  0,  (  «  rather  than  at  o  ■  0, 

^  ■  1  (which  is  no  longer  even  a  singular  point  of  the  differential 
equation).  Along  the  Integral  curve  passing  through  the  shock,  the 
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varlable  5  decreaBcs  from  Its  value  of  at  the  shock  to  some  minimum 
value,  and  then  increases  hack  to  this  value  at  the  axis.  Since  the  mag¬ 
netic  field  behind  the  shock  Is  given  by  Eq.  (24),  the  location  of  the 
singular  point  implies  that  the  current  In  the  wire  has  Jumped  to 
times  its  initial  wlue,  a  solution  in  contradiction  to  the  requirement 
of  a  constant  current  when  the  ambient  magnetic  pressure  is  retained.  The 
difference  in  location  of  the  singular  point  also  leads  to  discrepancies 
In  the  behavior  of  the  other  physlced.  variables  near  the  axis.  In  Pal's 
solution,  the  density  and  velocity  go  to  zero  at  the  axis,  but  with  a 
dependence  on  t|  which  is  qxiite  different  from  that  given  by  Eqs.  (27a)  and 
(27b).  The  most  dramatic  difference,  however,  is  in  the  pressure  which, 
rather  than  becoming  infinite,  as  given  by  Eq.  (27c),  goes  to  zero  at  the 
axlsl  The  solution  in  the  neighborhood  of  the  axis  thus  shows  a  sensi¬ 
tivity  to  the  shock  conditions  which  is  more  than  a  little  surprising. 


-20- 


VI.  MODIFIED  HYDRODYNAMIC  CASE 


The  modified  hydrodynamic  case  has  been  considered  in  detail  by 
Greenspan,  who  obtained  the  solution  by  the  numerical  integration  of 
three  coupled  differential  equations  combined  with  one  integral  of  the 
motion.^  The  solution  can  also  be  obtained  within  the  framework  of  the 
present  formalism.  Some  properties  of  the  solution  of  Interest  for  com¬ 
parison  with  the  pure  MHD  case  and  for  the  energy  considerations  to  follow, 
are  briefly  noted  at  this  point.  In  this  case,  the  magnetic  field  is 
assumed  to  be  continuous  across  the  shock;  the  boundary  conditions  at  the 
shock  are  therefore  simply  the  ordinary  hydrodynamic  strong-shock  condi¬ 
tions,  and  the  effect  of  the  magnetic  field  appears  only  in  the  flow  behind 
the  shock.  The  range  of  the  Alfv<^n  number,  is  now  Indeed  unbounded, 
in  contrast  with  the  pure  MHD  case. 

When  the  strong  shock  conditions,  together  with  the  continuity  of  the 
magnetic  field,  are  applied  to  Eqs.  (l6)  and  (17)»  one  obtains 


b 


2  ®o 


(28a) 


K ,  toJIr b 

(r-1) 


(28b) 


(ef.,  Eq.  (23)  tor  the  pure  MirD  case).  The  singular  point  of  Eq.  (l8) 
which  corresponds  to  the  axis  is  the  point  o  ■  0,  5  ■  •  Along  the 

integral  curve  through  the  shock,  the  parameter  ?  decreases  from  its  value 
at  the  shock,  ,  to  some  minimum  value,  and  then  increases  back  to 

^ThlS  was  only  a  special  case  of  the  more  general  problem  of  finite 
conductivity,  treated  by  Greenspan,  where  the  entropy  and  frozen  field 
Integrals  do  not  exist. 
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thia  value  at  the  axis.  In  the  modified  hydrodynamic  case,  the  magnetic 
field  behind  the  shock  is  given  by 

(cf.,  Eq.  {2k)).  Thus,  the  field  is  everywhere  smaller  between  the  shock 
and  the  axis  than  it  was  initially;  the  field  is  expanded  by  the  shock 
rather  than  compressed  by  it  as  In  the  pure  MHD  case.  At  the  axis,  the 
field  again  maintains  its  initial  value,  so  that,  as  In  the  pure  MHD  case, 

a  requirement  for  similarity  is  that  the  external  circuit  maintain  a  con- 

«• 

stant  current. 

This  difference  in  the  variation  of  the  magnetic  field  implies  a 
marked  difference  between  the  two  cases  in  the  current  distributions  in 
the  shocked  gas.  In  both  cases,  the  field  ahead  of  the  shock  is  assumed 
to  remain  unaltered;  the  integrated  current  through  the  gas  must  therefore 
be  zero  in  both  cases.  In  the  pure  MHD  case,  where  the  field  Jumps  across 
the  shock,  there  is  a  cylindrical  current  sheet  along  the  shock  front, 
flowing  in  the  same  direction  as  the  axial  current,  and  balanced  by  opposite 
currents  in  the  remainder  of  the  gas.  In  the  modified  hydrodynamic  case, 
where  the  field  Is  continuous  across  the  shock,  there  is  no  current  sheet. 
Between  the  shock  and  the  point  where  the  variable  (  reaches  its  minimum 
value,  there  are  currents  in  the  gas  flowing  in  the  same  direction  as  the 
axial  current:  these  are  balanced  by  the  opposite  currents  in  the  remainder 
of  the  gas. 

Although  the  singular  point  of  Eq.  (18)  for  Pal's  solution  is  the 
•axse  as  that  for  Greenspan's  solution,  a  current  Jump  is  implied  in  the 
former  case,  but  a  constant  current  in  the  latter.  The  reason  for  this 
im  that  pal's  Jump  condition  for  the  magnetic  field  leads  to  a  value  of 
b  given  by  Eq.  (25)  rather  than  Eq.  (26a),  and  hence  to  an  expression  for 
given  by  Eq.(2U)  rather  than  Eq.  (29). 
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Another  laportant  difference  between  the  two  cases  is  in  the  behavior 
of  the  other  physical  variables  In  the  neighborhood  of  the  axis.  The 
density  and  velocity  again  becone  zero  at  the  axis,  but  not  with  the 
logarlthale  dependence  on  t]  of  Eqs.  (27a)  and  (27b).  The  pressure, 
however,  also  goes  to  zero  at  the  axle.  In  contrast  to  the  pure  MHD  case 
where  It  becomes  Infinite.  The  results  for  these  two  cases,  as  well  as 
for  the  ordinary  blast  wave,  are  susmarlzed  In  Table  1. 
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Table  1 

DEPENDENCE  OF  PHYSICAL  VARIABLES  ON  NON-DIMENSIONAL  VARIABLE  r\ 
-  — ^IN'THE  NEIGHBORHOOD  OP  THE  AXIS  (n  -  0) 


VII .  ENERGY  COmCT 


The  one  remaining  parameter  of  the  problem  is  the  explosive  energy 
per  unit  length,  E^,  released  along  the  axis  at  t  =  0.  This  parameter 
enters  through  the  fact  that,  vlth  the  assumed  similarity,  the  total 
energy  between  any  two  similarity  lines,  in  particular  between  the  shock 
and  the  axis,  is  constant  in  time.  In  the  ordinary  blast  wave  solution^ 
the  initial  energy  of  the  gas  is  assumed  negligible  and  the  gas  is  assumed 
isolated  from  all  external  energy  sources.  The  energy  content  (entirely 
mechanical,  in  this  case)  of  the  gas  in  motion  is  therefore  equal  to  the 
energy  released  at  t  ■  0.  This  equality  allows  one  to  obtain  a  relation 
between  the  shock  speed  and  the  quantity  E^. 

When  a  magnetic  field  is  present,  the  situation  is  not  quite  so 
simple.  One  minor  point  is  that,  although  the  mechanical  energy  of  the 
gas  is  initially  negligible,  the  magnetic  energy  is  not  (being  infinite, 
in  fact)  and  must  therefore  be  taken  into  account.  Secondly,  while  the 
total  mechanical  energy  of  the  gas  in  motion  is  still  constant  in  time 
(as  is  also  the  magnetic  energy),  this  energy  can  no  longer  be  identified 
with  the  energy  released  at  t  *  0.  Because  of  the  coupling  between  gas 
and  magnetic  field,  some  fraction  (not  known  a  priori)  of  the  energy 
appears  as  mechanical  energy  and  some  fraction  as  magnetic  energy. 

Finally,  the  system  is  no  longer  Isolated  from  external  sources;  it  is 
coupled  to  an  external  circuit,  which,  moreover,  is  required  to  maintain 
a  constant  current.  It  is  conceivable  that  such  a  boundary  condition  on 
the  circuit  implies  an  exchange  of  energj'  at  t  •  0  between  external  source 
and  system.  That  this  is  In  fact  the  case  will  be  shown  below. 
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As  noted  above,  the  mechanical  energy  and  the  magnetic  energy  are 
both  constant  between  any  two  similarity  lines.  Another  quantity,  closely 
related  to  the  latter,  is  also  constant  between  any  two  lines  of  constant 
q,  viz.,  the  magnetic  flux.  Although  the  total  magnetic  flux  between 
shock  and  axis  Is^ constant  in  time  subsequent  to  the  explosion,  it  is  not, 
however,  equal  to  the  total  magnetic  flux  in  the  same  volume  prior  to  the 
explosion;  there  has  been  a  change  in  flux  (per  unit  length),  df,  in 
this  volume  given  by 


(30) 


Since  the  flux  between  shock  and  axis  is  constant  in  time,  it  follows  that 
an  amount  of  flux,  di,  is  supplied  (or  withdrawn)  at  the  axis  at  the  time 
t*0,  when  shock  and  axis  coincide.  Associated  with  this  flux  is  a  quan¬ 
tity  of  energy  (per  unit  length),  dE  -  I  dl,  which  is  supplied  to  (or 

c  o 

withdrawn  from)  the  system  at  t  ■  0  by  the  source  maintaining  the  constant 

current  1^.  The  increment  in  the  total  energy  of  the  system  thus  contains 

a  part  dE  in  addition  to  the  paz^  E  released  along  the  axis  at  t  ■  0. 
c  o 

The  increment  In  the  total  energy  of  the  system  appears  In  two  forms: 
an  increment,  dE^,  in  the  mecnanlcal  energy,  given  by 


(31a) 


*2  p  ^  r  p  1  d(6 

^^0  *  J  L  r^iT  "  2^ 

o 

and  an  increment,  fiE^,  in  the  magnetic  energy,  given  by 

4Eb  -  I  I  [t  ■ 

The  energy  balance  equation  for  the  system  then  becomes 

•E  -  -  AE  .  (32) 

0  0  B  c 

In  Eq.  (32),  the  quantity  AE^  -  AE^  represents  that  part  of  E^  vhich 
appears  as  magnetic  energy. 

The  energy  equation  can  be  written  in  the  form 

■  2no^c*^  InO„;r)  (5^) 

where  In(p^Jr)  is  a  constant  which  depends  on  the  parameter  and  the 
specific  heat  ratio  y,  and  which  can  be  evaluated  only  after  the  complete 

Integration  of  the  equations  of  motlor. It  is  this  equation  which  enables 

• . 

one  to  obtain  the  dependence  of  the  shock  speed  (which  is  related  to  c  ) 
on  the  explosive  energy  E^  and  the  Alfven  number 

^^In  his  treatment  of  the  modified  hydrodynamic  case,  Greenspan 
calculated  In(pQ;r)  only  in  the  limit  p^-  where  -  AE^. 
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VIII.  RESULTS  AND  COT?CLUSIONS 


The  results  of  the  numerical  calculation  of  the  energy  parameters 

are  shown  In  Flga.  5  and  6.  The  dependence  of  the  quantity  In(0^;7) 

on  Q  Is  shown,  for  two  different  values  of  7,  In  Fig.  5.  In  the  modl- 
*^0 

fled  hydrodynamic  case,  In(p^;7)  Increases  monotonlcally  from  thd  ordinary 
blast  wave  value  at  *  0  to  the  asymptotic  value  found  by  Greenspan 
at  pQ  ■*  ••  This  Implies,  from  Eqs.  (33)  and  (7),  that,  for  a  given  explo¬ 
sive  energy,  the' shock  speed  decreases  with  Increasing  magnetic  field.  The 
strength  of  the  shock,  however,  remains  constant  In  this  case.  In  the 

pure  MHD  case,  In(P  J7)  Increases  from  the  blast  wave  value  at  p  -0 

0  0 

to  a  maximum  at  a  0  of  about  0*2,  and  then  decreases  to  zero  at  3  ■ 

o  o 

For  a  given  energy  release  In  this  case,  then,  the  shock  becomes  weaker 

with  Increasing  magnetic  field,  but  the  shock  speed  first  decreases  to  a 

minimum  value  and  then  Increases  as  the  magnetic  field  Is  Increased. 

Fig.  6  shows  the  dependence  on  p^  of  the  fraction  of  which  appears 

as  mechanical  energy.  In  the  pure  MHD  case,  this  fraction  decreases 

monotonlcally  from  unity  at  p  «  0  to  l/r  at  p  »  1.  In  the  weak  shock 

0  '  o 

limit,  then,  there  Is  equipartitlon  of  the  explosive  energy  between  mechan¬ 
ical  and  magnetic  energy.  This  Is  In  contrast  to  the  modified  hydrodynamic 
case  where  this  fraction  Is  essentially  unity.  Independent  of  p^.  Actually, 
this  quantity  decreases  by  about  one  half  of  one  per  cent  at  a  value  of 
p^  of  about  0.3>  but  the  decrease  is  too  small  to  be  shown  in  the  figure. 

The  explosive  energy,  therefore,  appears  almost  entirely  as  mechanical 
energy  In  the  nodlfied  hydrodynamic  case. 

There  are  some  other  Interesting  differences  In  energy  content  In  the 
pure  VCTD  and  modified  hydrodynamic  cases  which  cay  be  Inferreu  from  Eqs. (30) 
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and  (51),  taken  together  with  Eqs.  (2li)  and  (29)  and  the  discussion  per¬ 
taining  thereto.  In  the  pure  MHD  case,  where  the  field  Is  everywhere 
compress'eS^behlnd  the  shock,  the  quantities  and  AE^  -  A*  are  both 
positive.  Both  quantities  are,  in  fact,  positively  infinite;  the  appear¬ 
ance  of  this  divergence  and  its  effect  on  the  validity  of  the  results  has 
been  discussed  in  the  Introduction.  The  difference,  AEg  -  AE^,  is,  however, 

finite  and  positive;  it  represents  the  portion  of  E  which  is  converted  to 

o 

magnetic  energy  at  t  ■  0.  Since  AEg  >  AE^,  this  part  of  E^  appears  in 
the  compressed  magnetic  field  between  shock  and  axis. 

In  the  modified  hydrodynamic  case,  where  the  field  Is  everywhere  ex¬ 
panded  behind  the  shock,  both  ABg  and  AE^  are  negative  (although  finite, 
in  this  ease).  The  difference,  AEg  -  AE^,  Is,  however,  again  positive, 
and  again  represents  the  portion  of  E^  which  is  converted  to  magnetic 
energy  at  t  ■  0.  However,  since  AEg  <  E^,  this  part  of  E^  does  not  appear 
now  in  the  magnetic  field  behind  the  shock  but  is  delivered  instead  to 
the  external  circuit. 

The  difference  Just  discussed  implies  a  corresponding  difference  in 
the  voltage  induced  in  the  external  circuit  by  the  gas  in  motion.  In 
the  pure  MHD  case,  AE^  is  positive;  i.e.,  the  external  source  supplies  a 
pulse  of  energy  at  t  ■  0.  Thus,  a  positive  voltage  pulse  (i.e.,  in  the 
same  sense  as  I^)  should  be  induced  in  the  circuit  at  t  ■  0.^^  Subsequently, 
as  the  shock  progresses  outward  through  the  gas,  ambf  .nt  magnetic  energy 
enters  the  region  behind  the  shock.  Since  the  total  magnetic  energy  (and 
magnetic  flux)  between  the  shock  and  the  axis  remains  eonstartln  time, 
magnetic  energy  must  leave  at  the  axis  at  the  sane  rate  it  is  entering  at 
the  shock.  Thus,  the  short  positive  voltage  pulse  should  be  followed  by 
a  long  negative  voltage  pulse. 
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In  the  modified  hydrodynamic  case,  the  situation  is  Just  reversed. 

The  quantity  dE  is  now  negative,  which  means  that  a  pulse  of  energy  is 
delivered  to  the  circuit  at  t  =  0.  This  implies  an  initial  negative 
voltage  pulse.  As  the  shock  progresses  outwards,  it  "pushes  out"  some 
of  the  magnetic  field  from  the  gas  through  which  it  passes  (the  field  is 
expanded  as  a  result  of  the  passage  of  the  shock).  There  is  thus  an  out¬ 
flow  of  magnetic  energy  through  the  shock  which  must  be  balanced  by  an 
Influx  of  energy  at  an  equal  rate  at  the  axis.  The  short  negative  voltage 
pulse  should  therefore  be  followed  by  a  long  positive  voltage  pulse. 

The  long  voltage  pulses  discussed  above  should  be  observable  experi¬ 
mentally.  However,,. whether  or  not  the  short  voltage  pulses  would  also  be 
seen  is  an  open  question.  As  mentioned  in  the  Introduction,  one  would 
expect  that,  in  any  actual  experiment,  where  the  axial  conductor  has  a 
finite  radius,  the  true  flow  would  approach  that  given  by  the  similarity 
solution  at  distances  from  the  axis  much  greater  than  the  radius  of  the 
conductor.  The  short  voltage  pulse  at  t  ■  0,  however,  is  a  property 
of  the  solution  which  is  associated  with  the  axis  itself,  and  one  may 
reasonably  wonder  whether  this  feature  of  -'he  solution  would  be  retained 
in  an  actual  experiment. 

Flns.lly,  the  energy  and  flux  considerations  outlined  above  help  to 
account,  at  least  heuristically,  for  the  difference  in  the  pressure  distri¬ 
bution  at  the  axis  in  the  two  cases.  One  can  think  of  an  imaginary  inter¬ 
face  at  the  axis  across  which  there  must  be  a  pressure  balance  if  the  flow 
la  not  to  separate  from  the  axis.  In  the  pure  mHD  case,  an  infinite 
magnetic  pressure  is  produced  on  one  side  of  this  Interface  by  the  addition 
of  magnetic  flux  at  t  ■  0;  this  is  balanced  by  an  infinite  gas  pressure 
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at  the  axis.  In  the  modified  hydrodynamic  case,  a  "partial  magnetic 
vacuum"  Is  created  at  the  axis  by  the  removal  of  flux  at  t  >  0;  pressure 
balance  In  this  case  requires  that  the  gas  pressure  be  zero  at  the  axis. 


sloral  variable  q,  for  various  valdes 
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Flg*  5  D«p«ndence  of  the  shock  speed  pemoeter,  In(fl  iy)t 
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